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Abstract 

To each 4x4 matrix of reals another 4x4 matrix is constructed, the so-called associate matrix. This associate 
matrix is shown to have rank 1 and norm 1 (considered as a 16D vector) if and only if the original matrix is 
a 4D rotation matrix. 

This rank-1 matrix is the dyadic product of a pair of 4D unit vectors, which are determined as a pair up to 
their signs. 

The leftmost factor (the column vector) consists of the components of the left quaternion and represents 
the left-isoclinic part of the 4D rotation. The rightmost factor (the row vector) likewise represents the right 
' quaternion and the right-isoclinic part of the 4D rotation. 

, Finally the intrinsic geometrical meaning of this matrix-based proof is established by means of the usual 

similarity transformations. 



y—i Notations 
> 

On Bold-face symbols: Geometrical entities 

^nI^ Normal-face symbols: Their algebraic representations 

' A: Van Elfrinkhof matrix (defined in Paragraph 4.2); Arbitrary 4D rotation matrix 

A/: Associate of matrix A (defined in Paragraph 4.2) 

' Mi, Mr: Matrices representing left- and right-multiplication by a unit quaternion, respectively 

(introduced in Paragraph 4.1) 

L = a + bi + cj + dk: Unit quaternion appearing as a leftmost factor in quaternion multiplication 

P = u + xi + yj + zk: Arbitrary 4D point represented as a quaternion 

R — p + qi + rj + sk: Unit quaternion appearing as a rightmost factor in quaternion multiplication 

5, Sl, Sr: Similarity transformation matrix and its left- and right-isoclinic components 

■ (Similarity transformations restricted to rotations of coordinate system) 

, (introduced in Section 5) 



1 Introduction 

The classical quaternion representation theorem for rotations in 4D Euclidean space states that an arbitrary 4D 
rotation matrix is the product of a matrix representing left-multiplication by a unit quaternion and a matrix 
representing right-multiplication by a unit quaternion. This decomposition is unique up to sign of the pair of 
component matrices. 

In this paper this theorem is proved as a result in the theory of matrices by pure matrix means. As a matrix- 
theoretic theorem it has no a priori geometrical meaning; for this one has to study the behaviour of its matrix 
formulation under a predetermined class of similarity transformations. 

The proof in this paper is not the first one in the existing literature; presumably Bouman ([BOUM 1932]) was 
the first person to publish a satisfactory proof of the representation theorem. 
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The application of quaternions to geometry is a well-known classical subject, but in the opinion of the au- 
thor it is also still a fertile research area. Think only of the use of low-dimensional geometry and topology in 
certain areas of theoretical physics. Or think of the hidden 4D rotational symmetry of the Coulomb field. 
The author's desire to make computer software for 3D and 4D geometry led to the purely matrix-based proof of 
the quaternion representation theorem for 4D rotations presented in this paper. 

2 Theorem and outline of proof 

Theorem: Each 4D rotation matrix can be decomposed in two ways into a matrix representing left -multiplication 
by a unit quaternion and a matrix representing right-multiplication by a unit quaternion. These decompositions 
differ only in the signs of the component matrices. 

Outline of proof: Let Af^, Mj^ be matrices representing left- and right-multiplication by a unit quater- 
nion, respectively. Then their product A — AIi^Mj^ is a 4D rotation matrix. 

Matrix Af^ is determined by four reals a, b, c, d satisfying the relation + 1)^ +c? +<P = 1. Likewise, matrix 
is determined by four reals p, q, r, s satisfying the relation -\- + r'^ + = 1 . 

The 16 products ap, aq, ar, as, . . . , dp, dq, dr, ds are arranged into a matrix AI , which has rank 1 and is easily 
expressed in the elements of A. Let us in this paper denote it as the associate matrix of A. 

Conversely, given an arbitrary 4D rotation matrix A, one calculates its associate matrix M in the hope that it 
is a matrix of products ap, aq, ar, as, . . . , dp, dq, dr, ds which are not all zero. This hope is vindicated by 
proving that M has rank 1 whenever A is a 4D rotation matrix. 

The proof is completed by observing that the sum of the squares of the elements of AI is unity, and concluding 
that two pairs of quadruples of reals a, b, c, d; p, q, r, s exist satisfying + 6^ + + = 1 , p'^ + q'^ + r'^ + s'^ = 1 
and differing only in sign. 

3 Quaternions in 4D Euclidean geometry 

In this paper the 4D Euclidean space is provided with a Cartesian coordinate system OUXYZ. Points are 
represented as column vectors {u,x,y, zY (in this paper denoted as the representation) or as quaternions 
u + xi -\- yj + zk. In this connection the coordinate system is also denoted as 01 UK. 

Both of these representations are based on an arbitrary choice of coordinate system. To obtain a possible ge- 
ometrical meaning of a matrix-algebraic result, one has to prove that it is invariant under coordinate system 
transformations {similarity transformations), which are in this paper restricted to rotations because we are inter- 
ested in Euclidean properties, not in more general afhne or projective properties. 

In the sequel we have to do with geometrical objects as well as their representations with respect to specific 
coordinate systems. For easy reading object and representation are denoted by the same symbol; the object itself 
in bold face, its algebraic representations in normal face, with decorations added as needed. 

4 Proof of the representation theorem for 4D rotations 

Refer to Section || above for the statement of the theorem and the outline of its proof. 
4.1 Isoclinic 4D rotations 

Let P be an arbitrary 4D point, represented as a quaternion P = u + xi + yj + zk. Let L = a + bi + cj + dk and 
R = p + qi + rj + sk he unit quaternions (a^ + b"^ + + d"^ = I, p'^ + q'^ + r'^ + = 1). 

Consider the left- and right-multiplication mappings : P LP and Mu : P PR. In the i?"' representation 
Ml and Mji are linear mappings with matrices 
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One easily proves that both Ml and Mfi are orientation-preserving isometrics of i?^ with the origin of coordinates 
O as a fixed point, i. e. rotations of R'^ about O. 
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Both Ml and Mn have the property of rotating aU half-hnes originating from O through the same angle (arccos a 
for Ml and arccosp for Ma); such rotations are denoted as isoclinic. There exists however a subtle difference 
between Ml and AI^, which is best illustrated by a specific example: 

Let Ml and M^ be left- and right-multiplication by the quaternion Qa = cosa + isina, respectively. Then Ml 
acts in both coordinate planes II and JK as a rotation through the angle a, while Mr acts in the II plane as a 
rotation through a and in the JK plane as a rotation through ~a. 

The two kinds of isoclinic rotations are conveniently distinguished as left- and right-is oclinic. 

Conversely, an isoclinic 4D rotation about O different from the non-rotation / and from the central reversion — / 
is represented by either a left-multiplication or a right -multiplication by a unique unit quaternion and so is either 
a left- or a right-isoclinic rotation. This theorem is not used in the sequel, but it allows us to speak of "left- or 
right-isoclinic" instead of the cumbersome "represented by a unit quaternion left- or right-multiplication". 
This theorem is presumably due to Robert S. Ball; in [BALL 1889] the author does not mention it explicitly as 
a theorem, but nevertheless gives a proof. However, Ball's proof covers only the case in which the off-diagonal 
elements of the given matrix are all nonzero. A complete proof is given in [MEBI 1994]. 

In Section ^ we prove that the left- and right-isocliny properties and the rotation angle are independent of 
the choice of coordinate system. 



4.2 General 4D rotations 

Let us study the composition of a left-isoclinic and a right-isoclinic rotation. First apply Ml, then Mu to a 
4D point P: one obtains {LP)R in quaternion representation. First apply Mr, then Ml'- one obtains L{PK). 
Quaternion multiplication is associative. Therefore left- and right-isoclinic rotations are commutative, and we 
have MlMrP = M^MlP representation) = LPR (quaternion representation). In the i?^ representation the 
mapping P LPR has the matrix 



MlMr = 



ap — bq — cr - 
bp + aq — dr 
cp + dq + ar 
dp — cq + br - 



ds —aq — bp + cs — dr 

- cs —bq + ap + ds + cr 

bs —cq + dp — as — br 

as —dq — cp — bs + ar 



—ar — bs — cp + dq —as + br — cq — dp 

—br + as — dp — cq —bs — ar — dq + cp 

—cr + ds + ap + bq —cs — dr + aq — bp 

—dr — cs + bp — aq —ds + cr + bq + ap 



(2) 



Van Elfrinkhof ([ELFR 1897]) was apparently the first person to treat the relation between quaternion 
multiplications and 4D rotations in this algebraic way. For this reason a matrix of the form of Eq. ^ with 
a'^ + b^ + c^ + d^ = 1 and p'^ + q'^ + r'^ + = 1 is in this paper denoted as a Van Elfrinkhof matrix. 



A Van Elfrinkhof matrix 



A = 



floo floi ao2 ao3 

flio ail 0-12 113 

020 121 ^22 123 

030 ^31 032 ^33 



(3) 



is readily interpreted as a set of 16 linear equations in the 16 unknowns ap, aq, ar, as, . . . , dp, dq, dr, ds. With 
slightly more work than just an inspection of the plus and minus signs in Eq. |^ one obtains these unknowns, 
arranged as a matrix, which is in turn written as a dyadic product: 
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(5) 

Matrix M contains at least one nonzero element and consists of mutually proportional columns and therefore has 
rank 1. The sum of the squares of its elements is (a^ + b'^ + + (P){p^ + + + s'^) — 1. In this paper M is 
denoted as the associate matrix of the matrix A. 

4.3 The general rotation as a product of isoclinic rotations 

In this paragraph we prove that any 4D rotation matrix can be decomposed into a left- and a right-isoclinic 
rotation matrix in two ways, differing only in sign. If we can find unit quaternions L = a + bi + cj + dk and 
R = p + qi + rj + sk such that the given 4D rotation matrix can be written as a Van Elfrinkhof matrix, then we 
are done. 

Let A (Eq. ^ be an arbitrary 4D rotation matrix, and let M, calculated according to Eq. || be its associate 
matrix. 

Easy and somewhat tedious calculations show that the sum of the squares of the elements of M is unity, and 
that all 2-by-2 minors of M are zero. In these calculations one uses the orthogonality of A and the equality of 
complementary 2-by-2 minors of 4D rotation matrices. 

The equality of complementary 2-by-2 minors follows from the general theorem stating that complementary 
minors of an orthogonal matrix B of arbitrary order are equal or opposite according to Det B equalling +1 or 
-1. ([CESA 1904], Art. 71, p. 54) 

It follows that M has rank 1, and there exist reals a, b, c, d; p, q, r, s satisfying + b'^ + + d^ — 1, 

+ + + = 1 such that M actually obeys Eq. § 
This yields two pairs of unit quaternions L — a + bi + cj + dk; R = p + qi + rj + sk, differing only in sign. 
Finally, Eq. gives the left- and right-isoclinic components, and Eq. || the original 4D rotation matrix as a Van 
Elfrinkhof matrix. 

Remark: It is sufficient to prove that the second, third and fourth columns of M are proportional to its first 
column. So only nine 2-by-2 minors need to be calculated, not all 36 of them. The Appendix shows one of these 
calculations. 

5 Independence of choice of coordinate system 

In this section we prove that the properties of left- and right-isocliny and the rotation angle of an isoclinic matrix 
are invariant under arbitrary rotations of coordinate system. This section is in part a refresher on dispacement 
transformations, similarity transformations, and their interplay. 

Let A be a 4D rotation with centre O. Let A be its matrix referred to a Cartesian coordinate system OUXYZ. 
Let OU' X'Y' Z' a be second Cartesian coordinate system, originating from the first one by a rotation S about 
O. Rotation A is a displacement transformation; rotation S is a similarity transformation. Let S be the matrix 
of S referred to OUXYZ and let Sl, Sr be its left- and right-isoclinic components. It does not matter which 
one of the two possible decompositions is taken. 

Furthermore, let p be an arbitrary 4D point, p — (u, x, y, z)'^ , p' — (u', x' , y', z')'^ its coordinates referred 
to OUXYZ, OU' X'Y' Z' , respectively. Then p' = S^^p expresses the new coordinates of p in the old ones, and 
p = Sp' the old coordinates in the new ones. Let Q = Ap point p rotated by A. 

How does similarity transformation S affect rotation matrix Al To find this out, begin at new coordinates p' , 
transform to old coordinates p = Sp', apply rotation matrix A to obtain Q — Ap — ASp', and finally transform 
back to new coordinates, ending at Q' = A'p' = S~^ASp'. We conclude that the matrix A' of A referred to 
OU'X'Y'Z' equals A' = S'^AS. 

Now we come to prove the intrinsic geometrical nature of the isocliny properties of 4D rotations. Consider 
the general rotational similarity transformation of the general 4D rotation: 
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A' — S ^AS — Sr ^Sl ^ -AlAji.SlSb^ — Sl ^A^Sl-Sr ^ArSr, 



(6) 



where the commutativity of left- and right-isocUnic matrices was apphed to sweep all left-isoclinic components 
to the one side while retaining their order, and the right-isoclinic components to the other side while retaining 
their order too. 

Applying Eq. |^ to a left-isoclinic matrix Al (set Ar = /) one obtains the transformation formula 

A'i^^S-^AlS = Sl-'AlSl, 

which shows that A^, being the product of left-isoclinic factors, is left-isoclinic. One may attribute the property of 
left-isocliny to the rotation Aj^ represented by matrix Al and all its transforms. This establishes the geometrical 
nature of this property. In other words, one may speak of left-isoclinic rotations in their own right, irrespective 
of particular matrix representations. 
Likewise, for a right-isoclinic matrix one obtains 

A'r = S ^ArS — Sii ^ArSb, 

which shows that A'j^ is right-isoclinic. Here too the property of right- isocliny may be attributed to the rotation 
Aj^ represented by Ar and all its transforms. 

In any rotation matrix of any order the main diagonal elements are the cosines of the angles through which the 
coordinate axes are rotated. Therefore, in an isoclinic matrix of either kind the main diagonal elements are equal. 
Their sum, the trace of the matrix, is invariant under similarity transformations. It follows that the rotation 
angle of an isoclinic matrix is also an invariant. So we may speak of the rotation angle of an isoclinic rotation, 
irrespective of particular matrix representations. 
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Appendix 

The sum of the squares of the elements of M in Eq. || is found to be 

4(aoo + «oi + «02 + ao3 + <4o + «!! + ai2 + ai3 + 020 + "21 + ^22 + + «30 + «31 + ^32 + a33)/16, 
which equals unity because of the orthogonality of A. 

All 2-by-2 minors of M are readily shown to be zero. For instance, 

16(mooTOii - miomoi ) = (aoo + an + + ass) (-ooo - an + 022 + 033) - {aw ~ oqi + 032 - 023) (oio - aoi - 032 + 023) 

= (022 + a33)^ - (aoo + flii)^ - (oio - Ooi)^ + (a32 - a23)^ 
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= + + ^32 + 13.23 ~ ^-Qo ^ '^11 ^ '^10 ^ '^Qi + 2(a22a33 — 0,320.23 — aooiii + aioOoi); (7) 

the bilinear terms cancel because in the rotation matrix A the complementary minors a22a33 — a32a23 and 
aoooii ^ oioaoi are equal. 

As regards the quadratic terms, we have 




therefore the quadratic terms in Eq. M cancel too. 
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